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Abstract 

We extend the Lagrangian formulation of relativistic non-abelian fluids in group 
theory language. We propose a Mathisson-Papapetrou equation for spinning fluids 
in terms of the reduction limit of de Sitter group. The equation we find correctly 
boils down to the one for non-spinning fluids. We study the application of our re- 
sults for an FRW cosmological background for fluids with no vorticity and for dusts 
in the vicinity of a Kerr black hole. We also explore two alternative approaches 
based on a group theoretical formulation of particles dynamics. 
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A de Sitter Group Generators [l8 
B Keplerian Rotational behavior \l9 



1 Introduction 

Recently there has been an enormous activity on fluid dynamics especially after the 
discovery of the Quark-Gluon Plasma (QGP) at RHIC. QGP describes a fluid at strong 
coupling regime of QCD where the usual perturbative description of fleld theories fail. 
There is now considerable evidence that it has several properties of a fluid and other in- 
teresting properties like chiral magnetic effect etc. For these reasons a self contained non- 
perturbative description of fluids is necessary. One of such attempt is due to AdS/CFT 
correspondence which naturally describes strongly coupled field theory. An advancement 
along this line from the perspective of Fluid/Gravity correspondence was developed by 
the work of Bhattachryya et al in pQ. In this article we discuss relativistic fluid dy- 
namics from the alternate view point of group theory. This work was developed in [2]. 
For a recent review see [3]. This technique relies on symmetry principles which makes 
it efficient to answer questions fike anomalies etc. as recently discussed in [4]. Such 
formalisms to describe particle motion were known from a long time (for some earlier 
references see [5]). The motivation for this formulation is not hard to understand either. 
To describe dynamics of particles or fluids in a general spacetime, it's natural to look 
at the symmetries that (quantum) gravity offers. Though the natural selection is the 
Poincare group (section l3.ll) . the incorporation of cosmological constant formally calls 
for the de-Sitter group in 4 dimensions. The formalism in one word gives a way to couple 
the particle/fluid with these gauge groups and deliver interesting results about them. In 
this formalism it also becomes easier to introduce fluid vorticity. 

The plan of the paper is as follows. For the rest of this introduction we discuss some 
models to describe relativistic abelian spinless fluids in the group theoretical language 
as was briefly addressed in [1]. In section [2] we describe the main model for describing 
relativistic spinning fluids based on a reduction of the de Sitter group. Later we present 
some alternate ways to formulate the problem and discuss each of their advantages and 
disadvantages. In section [3l2] we use a similar procedure using dS group, but in the con- 
text of torsion. Although torsion couples naturally with the spin density it is a property 
of spacetime that has not been observed; nevertheless this case should not be ignored. 
The model based on the Poincare group described in section 13.11 shows drawbacks over 
the others which come in the form of extra constraints that are not present in the other 
models. The fact that we can have several models is the generality of this approach. We 
use the formalism in situations interesting for astrophysical observations. For example 
in section [273] we describe the fluids in an expanding FRW spacetime. On the other hand 
section 12.41 is devoted to the study of dust in the spacetime of a Kerr black hole. We 
can use our formalism towards understanding some observable signatures for rotating 
galaxies and in particular in appendix |B] we describe their relation with Keplerian orbits. 
Appendix |A] collects some well known properties of dS group generators according to 
our convention. 
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1.1 Ordinary Abelian Fluid 

In 3+1 dimensions (with a mostly minus metric signature), an ordinary abelian fluid 
is described (in terms of collective variables and not as many particle motion) by the 
Continuity Equation 



and Euler's eqauation 



aip(t,r) + V(pv(t,r)) = 



dt'v + (v ■ V)v = force 



where the force is the force acting on the unit volume of the fluid, p and v are the fluid 
density and the velocity. These classical equations can also be found using a Lagrangian 
formulation with a Lagrangian [2] 



C = -fa, + -pv^ - V 



(1) 



where = (p, pv) and 

a, = - {d,e + ad,(3) . (2) 

Here 6, a and /3 are single valued non-singular functions which vanish at infinity. In 
this parametrization V6 is the irrotational part of the fluid velocity field while V x a = 
Va X V/3 is its vorticity. The above Lagrangian ([1]) can also be written as 



S 



+ pa$ 



^pv^ + V 



(3) 



where v = V6' + q;V/3. This type of parametrization is called Clebsch Parametrization. 

The Clebsch parametrization can be achieved in a group theoretical languag^ using 
SU{1, 1) or SU{2% With an element of the group parametrized by 



1 


'1 U 




" 


9= r. ^ 

V 1 — MM 


U 1 




_ e-''l\ 



where u = u\i g ^ SU(1, 1) and u = —u if G SU(2), the Clebsch parametrization is 
produced as 



where 



uu 



a 



uu' 



/3 = -z In 



9 describes the compact direction of the U{1) subgroup generated by |(T3 and a and (3 
parametrize the SU{1, 1)/U{1) or SU(2)/U{1) coset. The action ([3]) can now be written 

as 



JiJ 
2p 



+ V 



^Ajuy group can work equally well in principle. The fact that SU {1,1) works in this case, can be 
thought of as due to the fact that it has right number of parameters. We thank A. Polychronakos for 
a discussion on this point. 

^If we quantize vorticity (Va x V/3), then we need to work with SU{2). 
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The relativist ic version of such an action is given by 

^ = -'1 [r^T {asg^%g) + Fin)] 

with = nu^ for a unit flow vector u^, so that j^j^ = n"^, where n is number density, 
and F{n) encodes the specific dynamics (equation of state). The four velocity is then 
proportional to a^. Note here that for free theory (no potential V) which describes dust, 
F{n) = n. 



1.2 Non-Abelian generalization 

The generalization of the entire formalism to non-abelian case (for simplicity described 
in terms of SU{2)) at one particle level is described by the Wong equation with action 

mm 

S = j Knii' + AIQ'^x.-iwTiia-ig-^g) 

where a indexifies the number of generators of the gauge group and i the number of 
components of a vector. Here = Ti^ga^.g'^t"'), = |cr". For many particles the 
generalization is easy. With A indexing the particles we have 



S = ~i j dt^ wxTr{a3gx ^gx) 



with different gx^s for different particles which is in line with its continuous generalization 
of g{x). In fact to generalize it to the continuous case, we have, A — x, — / d^^/v, 
wx/v — ^ p{x)- We then get 



S = -i j d'^xpTr^cr^g ^g). 



Taking this last term as the term responsible for the symplectic structure, we can write 
an action 

^ = y d'xfTiiasg-'D^g) - j F{n) + Sym 

where j'^ is the current due to the charged particle density in the fluid as described 
before and is again given by = (p, vp) = nu'^ with = 1. 
For a general gauge group G, the action is given by 

S = -^ j / F{nr,n2,...) + SYM{A) 

(«) 

Here -R'(j) are diagonal generators of G and = n^^^-^ with i indexing the rank of the 

group G. 



2 Spinning Fluids Using the de Sitter Group 

From the models reviewed in section [1] it would be natural to write down a similar action 
describing the dynamics of fluids in terms of the Poincare group where the information 
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about the spin would be contained in the Lorentz subgroup and the information about 
the mass current would be contained in the translational part. As we will see in section [5] 
this creates internal inconsistency solved by the introduction of a constraint that reduces 
the degrees of freedom of the dynamics. Moreover the theory of Gravity cannot be 
written in terms of a gauge theory of the Poincare group directly. Only starting from 
the de Sitter group and then using a reduction limit we can reproduce gravity as a gauge 
theory. With this idea in mind we propose a generalization of the spin fluid action 
introduced in |3] with the following action to describe the full dynamics of relativistic 
fluids with spin for a torsion free backgrounc^l: 

5 = - y det e zj['^)Tr(a(m)To5f~ V/,5f) + zjf^)Tr(a(,)T235(" V^^f) + F{n(^m),n(s)) (4) 

where g is a. S'0(4, l)-valued field and Q;(m) and a(s) are constants respectively for the 
mass and spin coupling. The choice of the Lie algebra elements of so(4, 1) is taken 
considering that we are interested in studying matter fluids (Tq generates a time-like 
velocity) with spin, generated by T23 (see appendix |X] for a brief description of the 
generators of so(4, 1)). 

Action (jlj) is invariant separately under g — )■ ge~'^'^°^° and g — ge~^'^'^^'^^^ with 70 and 
723 constants; the associated conserved currents are with 



^.J?) = (5) 



where i = {m, s}. 

The explicit form of the current is 



= -z^Tr[K(,)^?-V,(?] (6) 

where = dF/dn^i), K^^) = a^rn)To and = a(s)T23. 
The variation of the action with respect to g yields 

V, [j^^^gK^^^g^' + jl,^gK^s)g-') = (7) 
that states the conservation of the so(4, l)-current 

J' = + 3U3K(s)g-' (8) 

Using (E]) we can rewrite ([7]) as follows 

3U \9~'^,9. K^m)] + J(^,) \9~'V,g, K^^)] = (9) 

where we used (E]). 

The components (refer to appendix [X] for the notation for the normalized trace) of 
the gK^i^g-'^ 

Q^f = -2TT[gK^,)g-'T^^] 

are the non-abelian charges of the theory ([8]). These charges are the coefficients of 
Swnab in the variation of the action with respect to the frame fields. In this sense they 



^For a discussion about torsion free models, see [4]. 
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are interpreted as spin charges. In particular Q^^^ is the spin charge of the fluid, been 
generated by the spin degrees of freedom, and Q^*^) is related to the vorticity of the fluid, 
being related to the mass current of the fluid. We will discuss more about the relation 
between Q^^^ and the vorticity in a short while. 

Equations can be used to write a Mathisson-Papapetrou equation for fluids. We 
apply the covariant derivative and then we anti-symmetrize the two indices and contract 
one index with the current to get the following expression: 

JwV[. (^^('m) = -^Tr [j(^,)V. {9Ki^)9~') V,gg^'] - '-Tr [gK^.^g-'R.,] j^^) (10) 

Adding the two equations in ffTOj) . and using ([7]), we have the generalized Mathisson- 
Papapetrou equation for fluids 

1 % 
= ~l^^A9K{m)g'^Rup\jlra) " 2^A9K{s)g''^ Rufj^jls)- (H) 

Notice that only the Q^-^ components of the charge contributes to the right hand side. 

Equation ( ITT]) can also be derived from the conservation of the energy-momentum 
tensor of the fluid given below: 

T'"" = - [^M^M + - nnM,n^s))] g'" + + ^-^(VT^) 

+2V, (j[i)Qj:) + jT„)g?:)) + 2v, [j[;,)Qj:^ + jI^^qj^^) (12) 

Setting Q^j?*^ = (and = if we want to discuss only about the mass current), the 
energy- momentum tensor (fT2l) reduces to the one for an ideal fluid with 

P = F (13) 

P + p = P{m)n(rn) + P{s)n(s) = F(^rn)n(rn) + F^s)n(s) 

where the = = dU / dN(^i) are the chemical potentials with U = pV the internal 
energy and A''(j) the particle number. Notice that the Mathisson-Papapetrou equation 
( ITTll is also easily derived from the conservation of the energy-momentum tensor 

where we used the equality 

which is the consequence of the anti-symmetry of Q'^^ and R'^ [p^^^ = 0. 

In section [2?T] we will perform the reduction limit of these equations. This procedure 
will reproduce the equation of motion for the fluid. The study of the effect of higher 
order corrections and the consequences of the possible extension of the model to the full 
de Sitter group are currently in progress. 
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2.1 Reduction Limit 

For the reduction limit we are going to consider the following parametrization 

g = Ae"^^ = Ah with A G S0(3, 1) 



(14) 



that describes a subset of the de Sitter group. The Lorentz Group is reproduced in the 
i? — > oo limit. Also noticing that 



aO A -1 



gK^s)g-' ~ a^^s)AT'^A~' + o(^ 

we re-establish the right proportion replacing at^s) with ai^s)/R- 
The currents (El) reduce to 



O 



1 



(15) 
(16) 



4-F(m) R ° -^(m) 



Qa 

—a(m)TaoA~^V ^A 



m 



(s) 



Tr 



a 



is) 



L R 



T23A-^V„A 



O 



i?2 



(17) 
(18) 



So, we see that both the j(^)'s are of the order O (-^). Hence we get 

1 



rif 



R 



implying that also F ^ O (-^) and ~ O (-^) . On the other hand, the so(4, l)-valued 
current (IHI) reduces to 



A 



2R' 



rpO ^ "a rj-iaO 



A-i + .r.^AT^^A-i 



R 



O 



1 



R 



o 



R^ 



where v'^ = A"q. 

To evaluate the reduction limit of ([7]) we replace g with the parametrization (1141) and 
substitute it in iQ obtaining 

where we kept terms until the (9(l/i?^). Evaluating the trace we get the following 
equations: 



the (a ^ 0) component 



the T component 



(19) 



[e'{A~W,A\' - d,e'] - 2a(.)jf,)(A-V,A)°Vi = 0; (20) 
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• the component 

[9\A~'V,A)^ - d,9'] - 2«(.)j[:)(A-V,A)°V^ = 0; (21) 

• the component 

[e'{A~'V,A\' - d,e-'] - 2a(.)j|^)(A-V,A)°V^ = 0; (22) 

• the T^^ component 

2a(,4)(A-V^A)i3 = 0; (23) 

• the T^^ component 

-2a(.)jf,)(A-V,A)i2^0. (24) 
Notice that the T^^ and the T° component are zero at all orders having 



at 







Multiplying equations (pi]) . (121]) . and (12^ respectively by 6'i, 02, and 6'3, and adding 
them, we find 



3 



- 2"M^'m E ^m^"^-^ - 2awjf,)(A- V,A)°P^=^1 = 0. (25) 

a=l 

Moreover equation ( |T9]) can be rewritten as follows 

a(,„)jf^)(A-i)^e^[9,(e^A'^o) + T^^ef A^ = 0; 
the above relation suggests that 

V" = e^A^o (26) 

is related to the fiuid velocity field and for the case of dust with zero vorticity satisfying 
the relativistic Euler equation ( |29]) . it is precisely the fiuid velocity (as discussed in 
section 12.21) . 



From the definition we have seen that the charges Q"^) are proportional to the 
spatial part of 6a which appears in several of the above equations. Moreover, defining 
the vorticity as in [11] 



n 



(m) 



we fincn 



oJuf, = —^j(ra)fL] = --Tr[V^ {gKi^ra)g~^) ^ ^,99"^] - -Tr[5(A'(„)^~^i?^J = 

V"-(m) / 2 Z 

if Q^^) = 0. Hence we are going to consider the Q'^-^ along with the spatial part of O"- 
to be the "charges" of the vorticity of the fiuid. 



^Notice that this recovers the definition in terms of tlie Clebsch parametrization. 
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Finally the Mathisson-Papapetrou equation in the reduction limit reduces to 



Vr, 



J{m) 



-Tr 



2R 



(f). 
R 



(27) 



Notice that the equation is of order as the energy-momentum tensor. This 

implies that, before taking the limit we also need to require the gravitational side to be 
of order O (-^) . 

We are going to end this section by checking that, in the reduction limit, the two 
Casimirs of 50(4, 1) reproduce the correct conserved quantities. The first Casimir is 
given by 



C2 = -v.u^rij'^r 



— > 



4 4 

and is related to the conservation of mass for a particle, n(m) being the number density 
times the mass. The second Casimir is given by 

9 'Jfi Jp^abcdQ 'Jfi'a'b''Ju'c'£ — ^ ^ abcdUl {s)a'b'Vc'^ 

that is proportional to the module squared of the Pauli-Lubanski vector, confirming our 
interpretation of as the spin charge of the fluid. 



2.2 Dust limit 

From the energy-momentum tensor ( |T2l) . once we neglect the terms containing the vor- 
ticity and the spin-density, it is easy to recognize the density and the pressure of an 
ideal fluid. 



that is. 



For spinless dust we can set 



(28) 



F{n(rn),n(^s)) = n^rn) and 0(5) = 0, 



that is, 



i(m)/. = -in(^m)TT[K(^rn)g ^V^^] and 

The Mathisson-Papapetrou equation (127|) then reduces to 



0. 



J(m) 



^(m) 



-zTr 



(Qa \ 



This recovers the relativistic Euler equation for a fluid of dust with zero vorticity 

1 
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(29) 



in the case 9°" = {6 , 0, 0, 0). Therefore from 

Jim) = ^«{m)9^^0 = n(^m)U^ 

we can recover 6q once is known. Notice that such an expression is exactly the same 
as that of a pressureless fluid with zero vorticity in the Clebsch parametrization 

The remaining part of the group- valued field A can be obtained from (fT9l) and (129|1 . 
which yields (126|) 

Jim) -eaA 

This tells us that v°'e'^ = A^qC^^ is nothing else than the fluid velocity with respect to 
the local frame described by the frame field. Notice that the above statement is always 
true in the case of dust with zero vorticity. 

Let us now consider the case of dust with spin. For dust with spin we need to set 

F(n(„,),r2(s)) = n^rn) + n(s)- 

The currents will be 



JM, = [d,9o + r (A-V,A),o] (30) 

Jis), = ^^(A-V,A),3. (31) 



nim) a(m) 

The Mathisson-Papapetrou equation (l27|l then reduces to 

^ ' "'^^M^M^'^o + ^W^^23^ A^: 



that recover the expected result (the Mathisson-Papapetrou equation for many non- 
interacting particles) in the zero-vorticity case, O"" = (^°, 0,0,0), 

2.3 Fluids in FRW spacetime 

In this section we are going to solve the fluid equations in an FRW background. Because 
of our initial choice of the action, only the matter dominated universe can be studied. 
The radiation dominated Universe would require a modified action with the correct 
requirement of producing a light-like current along with Tq — )> Tq ± Ti. 
We start with the general FRW-metric which is 



The static currents 



+r\d9^ + szn'9d(P^) 



1 — Kr^ 



^Mt = i^^[^t^o + r(A-^V,A).o] 
iRn^s) ais). I . 

J{s)t = -77^ VtA)23 

■iris) ri 
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depend only on the cosmological time for the case of an isotropic and homogeneous 
space-time under discussion. The conservation of these currents yields 







(32) 



with constant, and 



i(m) ( ^tQfm) + ~Qfm) ) + j{s) [ ^tQ(s) + 

a 



jU 9tQ'L + 2-Ql) +jU dtQZ + 2-Q 



a A«) 








where the Q(i)'s are only function of the cosmological time. 

The components of the energy-momentum tensor f|T2|) in our case are: 



= F + 4 



2 - kr' 



r 1 



4cot^(j*„,)g^„) 



et ^ 



2(2 - kr^ 



■^-U{rn)Q{m) +J{s)Q(s)) + ^(1 _ kr^^ U{m)Q"n) + i(.)<5( 



-2 cot ^(j(m)Q^m) + j\s)Q\s)j 



lls)Q%) 

We can also readily write Einstein's equations as 



a? k 



a 




Svr 



a 

a , 



F + 

2 



4(2 - kr^ 
r(l — kr"^^ 
k 

— - 87r(F(m)n(m) + F(^s)n(s) - F) 

2\ 



UUQlL) + + 4cot^(jf^)gf^) + j;,)Q 



a 



■t r^di 



(33) 



(34) 
(35) 



(36) 
(37) 

(38) 



The Ihs of equation (!36|) does depend on the cosmological time t only, while the rhs 
contains also other coordinates in the coefficient of the (5(i)'s. Equation (136 p is then well 
posed only if all the Q(i)'s, that are functions only of t, are zero or if 



Hm)Q^m) + Hs)Q' 



t f)tJ.V 



0. 



This second option is described by a proportionality relation between O'^A'^^^A'^^^ and 
A^pA^gj, that cannot be achieved. Therefore the only solution is to take all the (5(i)'s to 
be zero 



Q 
Q 



ab 

ab 
(m) 








a 



is) 







r = (^°, 0,0,0), 



that is, to consider a fluid with zero-vorticity and zero spin-density. This is not unex- 
pected considering that the FRW-background is homogeneous and isotropic. A more 
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interesting case will be to study inhomogeneous and/or anisotropic cosmological space- 
times. 

From the Einstein's equations we recover the classical FRW cosmologies. In partic- 
ular let us restrict to the case of dust where 



with Pq constant. We find: 



a? k ^ 
3-r + 3— = SttF 



a? a? 



fa k\ 



from which we can easily read ( l28ll 

p = F = poUi^m) and P = F[m)n(m) - F = 0. 

The constant po is just a scaling factor that can be reabsorbed by redefining n^s) and 
to set p = . In particular in equation ( l32i) we recognize the conservation of the 
rest mass during the cosmological evolution. 

We are now ready to find the group-fields in terms of 6^ and A. From the condition 
3{m) ~ ^(m) have that 

4i? 

— t + Co,o,o 



ml "5 



9" = 

L«{m) 

Considering that equation ( fT9l) reduces to 

we find that A"q = (1, 0, 0, 0) and therefore 

A = J ^ ^ with O e S0(3). 

This is not unexpected considering that taking a(^s) = enhance the symmetries of the 
initial action 

2.4 Fluids in Kerr-metric 

In this section we want to apply our fiuid description in a Kerr-background: 



, , , 2Mr\ , , 2Mr , ^ , , S , , ,.9/9 9 2Mr , „ , , 
ds^ = [1 — j dt^+2^-a sin^ 6dtd(l)-—dr^-me^- f + + sin^ 9 ) sin^ 



with 

A = + - 2Mr and S = + cos^ 6^. 

For the complexity of the equations we will consider the case of a stationary fiuid 
localized on the 9 = tx/2 plane 

f{i) = '^(i) Ni) (0,0,0, uf.)(r)] = n(i)M{i) (1,0,0, 
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and with zero vorticity 



Notice that such a current automatically satisfies the condition V^j|^^ = 0, n(j) being 
function of r only. 

Let's solve for dust. From the stationary condition and the expression for the mass 
current we have 



4i? 



that implies 
and hence 



= (^{m)t, 0, 0, 



9o{t,^) = Et-L(j) 



«(mV = ^(i?,0,0,-L). 



Here the constant E and L are respectively the energy and the angular momentum, per 
unit of n, of the fiuid (see appendix [B]) once we set Q;(m)/(4-R) to 1 with an appropriate 
choice of units. The above expression needs to satisfy the normalization condition 



fa(^m)E\^ lr^ + a^r + 2Ma{a-l).^ ^ ^ 



where we defined / = L/E, and where the angular velocity U(^m) = '^'(m)/'^\m) 

2Ma + (r - 2M)l 



r3 + aV + 2Ma{a - I) 



that in return implies that we are restricting to fiuids in rings around the rotating center. 
In appendix [B] we will explain how this indeed reproduces asymptotically the expected 
Keplerian rotational behavior in the galaxies. From now on, until the end of this section 
each expression is evaluated in r = rj, where rf is the radial coordinate of the particular 
ring under consideration. 

As in the previously studied case of dust with zero vorticity we can readily set 



\a _ a u 



or equivalent ly. 



AaO = ey(rn)tM 



AR 
4R 



e* -e^Z,0,0,e* -efO 

+ — a/ 1 , 

^= ,0,0, - a-0 

rVA r 



Using the equations (Ell), ([221) 

jf,)(A-i)^V^A'^o = with 6 = 2, 3 

reduces to 



(39) 



(40) 



(A-i)\ 



M .^a{M + r)\ 
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This equation has two possible solutions: we find uJ(^s) or we find A12 = A13 = 0. Here we 
study the second case that, as we will see, corresponds to not observe any spin current 
0). 

Using the group properties of the Lorentz group 



or.) 



A A" 







A -A". 



with i 



2,3 



and AaaA", 







we derive 



A, 



A 



30 



02 



A 



A32 
Ago 



Aoo 

A33 
Aoo 



32 



A - A 

^V03 — T ^i33- 

^^00 



1-(A 



A,5A, 



32J4-33 



23 J 



(A 



-AooA 



22^»-23- 



00 J 



= 1 - (A22)' 

From the above results we can finally write down 

/ Aqo A30 sin a A30 cos a \ 

_ 1 

COS a — sm a 

\ A30 Aoo sin a Aoo cos a j 

where Aoo and A30 are determined by equation (139|) . 

Notice that a rotation by a constant angle generated by T^^ redefines a. This is a 
symmetry of the action, so a is undetermined and can be set to zero to simplify the 
above expression 



A 



ah 



( Aoo 








A30 


\ 





1 
















1 







V A30 








Aoo 


/ 



Notice that the spin charges are not zero 

<523 oc a(s)Aoo 

while 

that satisfies the remaining constraint equations ( l23l) . ( 124|) . Having n(s) = simply 
means that the spin density ni^s)Q{s) is zero, that is, the distribution of matter in the 
Kerr-metric does not show any spin polarization. 



3 Alternative Approaches 

In this section we take up two alternative approaches to generalize the relativistic spin- 
ning fiuid action introduced in [1]. The first is based on a direct use of the Poincare 
group following [5]. The second approach uses again the de Sitter group and a reduction 
procedure similarly to section [2] but in the contest of Mac-Dowell-Mansouri gravity in 
which the coupling between spin and torsion comes naturally. Also this approach finds 
its origins in [3] but in a slightly different version described in [H] that makes it suited 
for Mac-Dowell-Mansouri gravity. 

One main difference between the models described in this section and the one in 
section [2] is that here both models contains only one single current while in the previous 
formulation we separate the spin and the mass transport into two currents, effectively 
increasing the degrees of freedom in the dynamics. 
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3.1 Fluids Using Poincare Group 

In this section, to find both the mass and spin transport by a fluid we use the full 
Poincare group. We generalize to fluids in line with the model introduced for relativistic 
spinning particles in j5]. 

Following |5], the action for spinning fluids introduced in |1] can be generalized to 



S = — det ed X 



(41) 



where u = (1,0,0,0), the function F depends on the equation of states and A has the 
units of angular velocity and represents the "spin-charge" of the particles composing the 
fluid. Notice that, like in the previously described models, is conserved, the action 
being invariant under the right action of a constant element in the Tu direction. The 
variation with respect to gives 



-F'{uA-')ael + z^Tr(Ti2A~ V^A) = 



or equivalently. 



where 



V, = ;^^Tr(T,2A-V,A) (42) 

V, = {uA~%el = (A"i)V;: 

is the fluid velocity fleld naturally constructed as boost with respect to the local frame 
fleld. The variation of the action with respect to A yields the spin-precession equation 

2F'j[^;^l - V^ij^S''^) = (43) 

where zAATi2A~^ = 5 is the spin-density. 

The Mathisson-Papapetrou equation can be constructed following the same steps as 
in section [2] 

f^uiF'v,) - fV,{F'v,) = ^Tr(j-V.5V^AA"i) + Ti{SR,,)f (44) 

where we multiplied by j''. 

The energy-momentum tensor is given by 



nap 



-F[{uA^\el3^] + z^j''Tr(Ti2A- V,.A) 



g'-P + F'vac'^^r - V^(j^g^" + j"Q^^) 
v^fF' - F) g"^ + F'yf^r - V^(j^g^" + rQ^^) 



Note that the energy tensor T"^ is not symmetric because of the appearance of the 
non-symmetric factor —F'v^j". 

The main reason for which in general we have to expect a non-symmetric energy- 
momentum tensor is the presence of a term in the action linear in the frame fleld: 
j^VaC'^^. Only in the case in which we have an action in terms of the frame fleld such 
that only quadratic combinations appear, we will have a symmetric energy-momentum 
tensor because each quadratic combination of the frame fleld can be replaced by the full 
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metric. The non-symmetry of the energy-momentum tensor imphe^ the presence of a 
constraint in the theory as is obvious from the Einstein's equation 



Indeed the Ihs of the above equation is symmetric by construction which in return imphes 

rp[»l3] ^ _p'^ll3ja] ^ (45) 

implying that and are proportional! 

The first consequence of the constraint f l45l) is that equation becomes 

V^ij^S'^") = (46) 

effectively limiting the dynamics. The second consequence comes straightforwardly con- 
sidering that the parallel-condition 

r = 



together with equation (j^2|) and the dispersion relation j°'ja = n^, imply 



2 2 a 2 

n = a VaV = a , 



that is, 
and 



^2 _ ,a 



This means that is the current of the fluid and that some of the degrees of A are 
constrained by the specific choice of F{n). Note that we want to get dynamics with F{n) 
in the action. To achieve it for mostly minus metric signature, we need to choose the 
plus sign in the above relation. Notice that the current was not previously defined 
by an equation of motion, its definition in fact comes from the above constraint. 

In the presence of our newly found constraint, the Mathisson-Papapetrou equation 
reduces to 

3'V,{F'v,)-3^V,{F'v,) = TT{SR,,)f 



3.2 Fluid Using de Sitter Group with Torsion 

Mac-Dowell-Mansouri gravity is a 5*0(4, l)-gauge theory with a term that breaks the 
symmetry down to the Lorentz group to recover Einstein's Gravity with torsion. In this 
framework it is natural to introduce particles within the group theoretical approach of 
[5] as described in |9]. We want to generalize this procedure to fluids. 

The fluid action will be of the same form of the actions that we described previously: 

S= - [dete [ifTT{Kg~^D^g) + F(n)], (47) 



second option would be a modification of the gravitational sector tliat would produce a non- 
symmetric Einstein tensor. In this paper we shall not address this option. 
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where g e 5'0(4, 1), K = cqTq + C23T23 and 



The gravitational potential above is the 5'0(4, l)-gauge potential of the Mac-Dowell- 
Mansouri Gravity where / denotes the de Sitter radius; / — 00 is the reduction limit of the 
theory. The variation of the full action, including the gravitational sector, with respect 
to the frame field e will produce Einstein's equations while the variation with respect 
to the spin connection w will give us the torsion equation. As we have already seen in 
the previous section, the presence of the frame field e has to do with the translational 
part of the Poincare; indeed in the reduction limit the coefficient Cq will correspond to 
the momentum density of the fluid. 

The gauge group is broken from the gravity sector to 5*0(3, 1) but this does not 
reduces the fleld g that will still have values in S'0(4, 1). 

Noticing that the right action of / = exp aT^ would send an element of the S'0(3, 1)- 
gauge subgroup into an element of the 5*0(4, l)-gauge group while / = exp/3T23 would 
leave it in the 5*0(3, l)-gauge subgroup and considering that an 5*0(3, l)-gauge trans- 
formation is internal in 5*0(3,1), we see that it is consistent to study the restriction 
of the (yf-fleld to the 5*0(3, l)-subgroup only. The general case of (7 G 50(4, 1) has a 
more rich dynamics but here we are going to restrict to the case 5*0(3, l)-valued g that 
reproduces the expected equations [9]. 

The conservation of the current follows from the variation of the action (H7|) with 
respect to g, giving 

= D,{fr') = d,{fr') + T>;^fr' + 2rw,y' - ^j^eiyi (48) 

where 

J = gKg-' = coKTa + casA^A^T,, = p'^Ta + s'^'Tab 

The action is also invariant under the right action g — )■ gf, with / G 50(3, 1) and 
commuting with K. It is easy to see that the only commuting elements in so(4, 1) with 
K are the generators Tq and T23. Such global symmetries imply the conservation of the 
current 

= 0. (49) 
Using ( H9|) . the previous expression (H8|) reduces to the spin precession equation 

r'^f.r' - ^J^ety] = (50) 
The variation with respect to the current produces 

J, = -^-^r[Kg~'D,g] 

This, following the usual procedure already deflned in section [2] give the Mathisson- 
Papapetrou equation for fluids 

zTt {fV,JV,gg-') - 2ifTi [JR,^) 
zi-Tr (j^eiy r,,V^(7^?-i) - 2y-Tr {JR,^) 
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4 Conclusion 



In this work, we give the Mathisson-Papapetrou equation for fluids looking to effects 
relevant in strong gravitational background^. Our approach rely completely on the 
group theoretical approach started with works of various groups, as mentioned in the 
introduction. In sections [2] and [3], we have presented several models based on different 
prescriptions reproducing the whole dynamics of a fluid with spin in terms of both the 
mass and spin transport. 

With the discovery of QGP in the laboratory and also for intergalactic fluids, such 
types of studies are relevant and especially in the case of QGP, the different types 
of models that we have presented could be put to direct tests. Moreover, using the 
AdS/CFT correspondences, we could also test these models using the well handled 
weakly coupled gravity description. 

As per future directions, there are several avenues to proceed. On the negative side 
of this formalism we note that it's quite difficult and we may have to incorporate the 
viscosity and other dissipative effects in an ad hoc way. But the formulation allows us to 
study various general backgrounds in the presence of fluids with spin densities making 
this interesting for astrophysical studies. Specially, as we saw, our formalism can easily 
account for fluids in a different spacetime. We chose to explore the Kerr black hole 
case as this might have an application in astrophysics, namely in the study of rotating 
galaxies and the rotating black hole in their centei0. 

In terms of the general formalism, based on the work of anomaly in [1] , one might 
try to explain the chiral magnetic effects in the RHIC experiments. Moreover this could 
be thought of as a parallel description of Fluid/Gravity correspondence or the work 
done by AdS/CFT community. Recently Bredberg et al. [10] have uncovered a possible 
connection between incompressible Navier- Stokes equations and a corresponding dual 
Einstein's equations in one higher dimensions. If correct, it would be interesting to see 
if that can be described from the formalism we have adopted here. 
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Appendix 
A de Sitter Group Generators 

In this appendix we list the convention used in this article for the representation of the 
Lia algebra associated with the de Sitter group. 

^The fact that it is in the realm of strong gravitational effect can be understood by the fact that the 
spin coefficient is small unless it is coupled to gravity strongly. Similar reasoning holds for Torsion. 
^Wc thank Matthew O'Dowd for teaching us some astrophysics at this point. 
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The generators of so(4, 1) are constructed, similarly to so(3,l), starting the four- 
dimensional Dirac's gamma matrices 

with the addition of 7^. Such algebra can be rewritten in a shortened notation as 

where A takes values {0, 1, 2, 3, 5} and 77"^^ = {1, —1, —1, —1, —1}. 

Similarly to the case of the Lorentz group the de Sitter group is generated from Lie 
algebra generator defined as follows 

Tab = ^ [7a, 7b] with T^s = -Tg^ = Ta 

o 

The commutation relations between the Lie Generators are 

2[T^_B, Tcd] = VbcTad — VacTbd — VbdTac + VadTbc, 

that is, 

2[Tafe, Tcd] = VbcTad — VacTbd — VbdTac + VadTbc 
2[Tab, Tc] = TjbcTa - VacTb 2[T„, Tj = Tac 

The normalized traces of the Lie generators are given by 

-2Tt[TabTcd] = ^{vacVbd - VadVbc)- 

B Keplerian Rotational behavior 

To simplify the problem we assumed from the beginning circular orbits. In this appendix 
we will also consider a non-zero radial component in the current and only at the end we 
will set such a component to zero to recover circular orbits. 
For the dust case under discussion we set 

j{m)^l = nu^ = n{E, Ur, 0, -L). 

The constraints on come from the zero vorticity condition and its normalization!^: 

V[^n^] = n^n^ = 1. 

Multiplying the zero vorticity condition by 2u^ and using the normalization condition 
we find 

= 2n^V[^n^] = uf^V^u^. 

Let's start from u^'Vfj,Ur- AH the other components will be automatically zero once we 
set Ur = 0. From the normalization condition 

1 = —EL —L^ - —ul 

rA rA rA 



^Wc have already shown that the current is trivially conserved for m,. = 0. 
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we have 



U^drUr = — - 



r 
A 



M 



rA 



Therefore we can easily expand 



1 / r 



2 V A 



4r' 



,r - M 



A2 



A2 



3r'^-4aV^+10Mr*-aV+6MaV^-8MV 



„4 „4 



2r2A2' 



2„2 o ;\/r2„3 



+ (3r^+4aV^-8Mr^+aV-6MaV2-2a^M+8M2aV)E2-(-r^+aV+6Mr2-8MV+2a2M)L2+ 

+4aM(3r^ - 4Mr^ + a^r)EL] 

In this appendix we are interested in showing that the theoretically expected Keple- 
rian rotational behaviour is reproduced, therefore we will consider only the asymptotic 
behaviour, for which we can take a — > 0: 



\ r . -M 
2 I A A2 



A' 2 



2r2A2 



-4Mr3(r-2M)+4r(r-3M)(r-2M)L2 



Setting Mr = to simplify our condition (in this way we are going to assume only circular 
orbits), we have 



Mr) 



and hence 



E 



rf-3M 



(r/ - 2M)2 
rfirf - 3M) 



that indeed reproduces the expected Keplerian behaviour of the rotation curve for galax- 
ies 



gt^E 




M 
J- 
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